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ABSTRACT 

Given  a  Jordan  curve  T  in  the  complex  plane,  we  describe  a  polynomial 
family  which  is  asymptotically  orthonormal  on  I*.  The  polynomials  have  some 
similarities  with  the  Faber  polynomials  but  are  simpler  to  compute  with. 
Numerical  examples  are  presented. 
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Work  Unit  Number  3  (Numerical  Analysis) 
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SIGNIFICANCE  AND  EXPLANATION 


Let  r  be  a  Jordan  curve  in  the  coaplex  plane.  Me  describe  a  polynomial 
family,  which  is  asymptotically  orthonoma  1  on  r  (as  the  degree  of  the 
polynoaials  increases ) ,  and  which  is  siaple  to  compute.  We  use  the 
polynoaials  to  approxiaate  functions  f(s)  analytic  in  the  interior  of  T 
and  continuous  on  I*.  After  soae  initial  calculations*  which  are  independent 
of  the  functions  to  be  approximated ,  each  functions  f(s)  can  be  approximated 
by  an  nth  degree  polynomial  in  0(n  log(n) )  operations.  Theoretically,  one 
could  use  the  Faber  polynoaials  for  T,  but  the  proposed  polynoaials  are  auch 
simpler  to  compute  with  end  almost  as  effective,  as  is  seen  by  the  bound  for 
the  resulting  polynomial  projection  given. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  NEC,  and  not  with  the  author  of  this  report. 


AM  ASYMPTOTICALLY  ORTHOMOAMAL  POLYNOMIAL  FAMILY 


Lothar  Reichel 


1*  IMTMOOOCTION 

Let  r  be  a  Jordan  curvo  tn  the  cooplox  plane,  and  let  a  -  4(w)  be 
function  on  |w|  >  1,  such  that  4  nape  |w|  >  1  onto  the  exterior  of 
4(«0  ■  •.  4  can  be  extended  to  a  continuous  bij active  step  on  |w|  >  1. 
we  will  study  are  defined  by 


n-1 


n-1 


(1.1)  P  (s)  I-  ~  (  n  (s  -  4(«aki’/n)>  +  B  (*  -  4<e<2*M)1*/,*)>), 

it  _  n  , _ _  , _ _  ' 

2c  k“0  k"0 


where  c  denotes  the  capacity  of  1*. 


an  analytic 
r,  with 
The  polynoaials 


n  •  0,1,2,..., 


■x.  1.1.  Let  r  be  the  unit  circle,  and  let  s  *  4(w)  -  w.  Then 

.  ui  - 1  (■;’  >. .  .  v -  .<*«»*.,, . 

k-0  k-0 

-  i  <(s*  ♦  1)  ♦  (s*  -  1)>  -  s",  n  -  0,1,2,..., 

These  are  the  Faber  polynosdals  for  the  unit  disk.  ■ 

Sx.  1.2.  Let  r  be  the  ellipee  {x  ♦  iy,  (*)2  ♦  (*)2  •  1 },  a  >  b.  Let  a  i«  and 

d  •*  (a2  -  b1)1^.  Then 

(1.2)  s  •  4(*>  ■  da(w  ♦  e'V1)  . 

The  capacity  of  T  is  da.  Substitute  (1.2),  wt  *■  ei,t/n  and  c  -  da  into  (1.1). 

This  gives  after  some  simplifications 

.  .  _  n  .  -2n  -n 

p„(s)  -  w  ♦  a  w  . 
n 

These  are  the  Faber  polynomials  for  the  ellipse  bounded  by  r,  c.f.  Curtiss  (2).  For 
future  reference,  we  note  that 

Sponsored  by  the  United  states  Army  under  Contract  Mo.  DAAQ29-80-C-0041 . 


(1.3) 


«  ♦  • 


PB(*)  ♦  w  , 

if  a  >  1,  or  equivalently,  if  b  >  0. 

Xn  Motion  2,  wo  show  that  for  a  large  claaa  of  boundariaa  r  tho  polynomials 
pn(s)  ara  asymptotically  orthonormal  with  r as pact  to  an  inner  product,  c.f.  (1.2)  and 
(1.3).  Wo  uaa  thasa  polynomials  to  dafina  a  bounded  projection  operator  for  polynomial 
approximation  of  functions,  which  are  analytic  interior  to  T  and  Mtisfy  certain 
smoothness  properties  on  F.  Xn  section  3,  we  show  how  this  projection  onto  polynomials  of 
decree  <  n  can  be  computed  in  0(n  log(n))  operations  for  each  function  to  be 
approximated,  provided  that  sosm  initial  calculations  independent  of  the  function  to 
approximate  have  been  carried  out.  Section  4  contains  numerical  axasqples. 


it  of  a. 


Let  0  denote  the  ooan  lntarior  of  r,  and  lat  0  ba  tha  ooaplwaen 

c 


Tbaoraai  3.1. 

amom  that 
|w|  -  1  for 


1 

dv 


is  continuous  on 
j  >  0.  Than 


M 


1. 


and  that 


is  of  bounded  variation  on 


(2.1)  p  (*)  “  wB(1  ♦  o( n”3 ) ) ,  n  ♦  uniformly  for  t  «  0  ,  vhara 

n  o 

w  is  defined  by  ♦(*)  »  a  . 

(2.2)  p_(s)  “  o(#*^*'),  n  ♦  •,  id,  and  uniformly  for  s 

n 

belonging  to  any  closed  subeet  of  0. 
If  +  is  analytic  on  |v|  -  1,  then  there  is  a  constant  r,  0  <  r  <  1,  such  that 


p  (s)  •  wB  ♦  0(r“B),  n  ♦  uniformly  for  s  «  0  ,  where 

n  o 

»  is  defined  by  w  -  ♦(«). 

p  (s)  “  0(r"*),  n  ♦  •,  stfi,  and  uniformly  for  s 

U 

belonging  to  any  closed  subeet  of  0. 

Proof,  in  the  proof  ms  make  use  of  results  Curtiss  [2]  obtained  in  his  investigation  of 
the  product 

V  i.  -  «.’•“'*»  . 

k-0 

Curtias  (21 >  Lanma  1 ,  shows  that  if  f(d)  is  a  2*~periodic  oamplex  valued  function  of  the 
real  variable  •,  absolutely  continuous  on  the  interval  0  <  6  <  2f,  and  if  ^  ■  2i 
k  “  0(1 )n,  then 


(2.3) 


__  n-1  2* 

»  I  «•„)  -  /  f( #)d8  ♦  ©(n"1), 

k-0  *  0 


d3f 


If  — r  is  absolutely  continuous,  Curtiss's  proof  of  (2.3)  supple 
dd3 


ated  by  integration  by 


parts  yields 


9  In’ 


~  l  f(8k)  -  /  f  ( 9)d9  ♦  o(n  3  ),  n  ♦  - 


Following  Curtins,  wo  introduce 


g(w,w)  »■ 


c(w  -  w) 

9*<w>  ,  w  -  w  , 


whoro  |w|  >  1,  |w|  -  1.  Lot  9(9, w)  «■  log(q(e  ,w)>.  With  a  branch  of  tho  logarithm 
choooon  so  that  w  ♦  9  is  analytic  and  single  valued  for  |w|  >  1,  continuous  on 
|w|  “  1  and  vanishes  as  |w|  ♦  ",  the  Cauchy  integral  fontula  yields 

21 

/  9(8,w)d9  ■  0,  |w|  >  1  . 

0 

Let  9fc  ««  2«k/n,  k  ■  0(1  )n  -  1.  Then,  with  s  -  9(w), 


n  (x  -  9(«  )> 

_o _ 

cB(wB  -  1) 


rl  nrl  18 

I  -  I  ,{v  \  . 


0  c<w  *  • 


With  a  suitable  branch  of  the  logarithms,  wo  have 


H  (s  -  9(e  ")) 


n-1  n-1 

r  m  .  _  n  r 


— )  •  £  "*•  V  •  b  £  «">  W.  -  v 

c  (w  -  1)  k*0  k»0 


-  i?  (I  ■  V  -  /  **■•>«•>  -  v«> 

k-0  0 


n-1  19 

—  9  (s  -  9(«  ))  -  (wn  -  1X1  ♦  0(H(s))),  n 

_n  .  -  n 

C  IC"0 


Curtiss  (1]  considers  the  case  J  -  0,  and  shows  that  R_(9(w>)  -  o(1),  new,  uniformly 

n 

for  |w|  >1,  For  j  >  0,  some  straightforward  modifications  of  Curtiss's  proof,  like 


replacing  (2.3)  by  (2.4).  yields  H  ($(*))  *  o(n”^),  n  ♦  •  uniformly  for  |w|  >  1.  If 

n 

4(«)  is  analytic  in  a  neighborhood  of  the  unit  circle,  then  (i  ia  an  analytic  function 
of  both  its  arguments,  and  |H  (+(w))|  <  Arn  for  some  constants  A  and  r,  0  <  r  <  1  as 

n 

n  ♦  •,  uniformly  for  |w|  >  1. 

How  replace  0  by  8  i-  0  ♦  k  -  0(1  )n  -  1,  in  (2.5).  Then 
K  JC  K  n 


(2.8) 


n-1 

8 

k*0 


i6k 

g(w,e  ) 


n-1 

n 

k»0 


16, 


c>) 


n.  n 
c  (w 


+  1) 


and  analogously  to  (2.7)  we  obtain 

1  *  lK  „ 

(2«»)  “T  8  (a  -  +(e  “))  -  (wn  ♦  1)(1  ♦  o(H  (a))),  n*». 

c*  k-0  B 


The  average  of  (2.7)  aad  (2.9)  yields  (2.1).  Also  (2.2)  follows  from  results  of 
Curtiss.  For  0  of  bounded  variation  on  |w|  »  1,  Curtiss  shows  that 


(2.10) 


1  *‘1  i#k 

~  *  (a  -  +(e  ))  -  -1  ♦  0(h(n) ) ,  h(n)  -  o(1),  n  ♦  •  , 

0  k-0 


for  any  a  8  fi.  and  uniformly  for  a  belonging  to  any  closed  subset  of  Q.  Again  it  is 

straightforward  to  show  that  if  — *  is  of  bounded  variation  for  some  j  >  0,  then 
-1-1  ** 

h(n)  -  o(n  ).  If  4  is  analytic  in  a  neighborhood  of  |w|  -  1,  there  are  constants 
B.  r,0<r<1.  such  that  |h(n)|  <  Br  n.  Inspection  of  Curtiss  proof  also  shows  that 


1  0-1  l,k 

<8.11)  “  J  (a  -  f(e  ))  -  1  ♦  0(h(n)),  n  ♦  •  , 

c  k-0 


for  any  a  «  0.  uniformly  for  z  belonging  to  any  closed  subset  of  Q.  Adding  (2.10)  and 
(2.11)  yields  (2.2).  e 


Ut  4  (*)  denote  the  Invars*  asp  of  ^(v).  Ha  introduce  tha  innar  product 


(2.12)  (f,g>  «-  ^  /  f(*)g(s)|d4_1(s)|  -  /  f(4(w))g(4(v))|dv|  . 

T  |w|-1 

Tha  bar  danotas  coaplax  conjugation. 

Thaoraa  2.2. 

(2.13)  (p  (s),**)  -  0,  i  >  n,  w  -  ♦_1(a)  . 

II 

(2.14)  (p  (*),*“)  -  1  . 

n 

Proof. 

(2.13)  /  V  (♦(*)-4(a2Wik/,1))wt|d*|  - /  V  ( ♦(*)-♦( a2 ’lk/B ))*"*" 1  da 

”  1*1-1  k-0  2n  |w|-1  k-0 


4  baing  analytic  axtarior  to  |w|  -  1,  w*  can  replace  tha  integration  path  by  a  circle 

|*|  -  R  sufficiently  larga  so  that  4  has  an  expansion 

(2.16)  4(*)  -  c*  aQ  +  a^*-1  ♦  ♦•••*  1*1  >  R  . 

Substituting  (2.16)  into  (2.13)  yields  (2.13).  (2.14)  follows  frcai 

k  1c-1  . 

Pk(4(«))  -  **  ♦  l  * 

J “  • 


Ha  will  use  a  Petrov-Galerkin  eethod  to  ooapute  polynomial  approxinations.  Let 
0<n>  denote  the  Sr  a— Ian 

“  lo^l  a  Cw  l-  (P4(s),*k),  0  <  k,t  <  n,  *  -  4-1(s)  , 

and  let  P*®*  be  the  Fourier  operator 

F(n)f  -  [F^J,  ?k  «-  (f(s),wk),  0  «  k  <  n,  *  -  4“1(s)  . 

TheoraM  2.3. 

bet  f(s)  be  analytic  in  8,  and  have  uniformly  bounded  Fourier  coefficients 
F.  -  (f(s),*k),  s  -  4(f)*  hot  the  projection  P_  be  defined  by 


with  a  »  a  «■  Him,  if  it  of  bounded  variation  on 

(ht 

Iwl  -  1, 


(2.1»)  ak  ♦  fk,  k  ♦  •»  *  >  k 

(a) 


froof.  by  (2.13),  (2.14)  0'"'  ia  an  uppar  triangular  matrix  with  diagonal  alaaants 

6^  •  1.  Tha  uppar  triangular  >1— anta  hava  by  thaoraw  2.1  tha  fora 
0tt  “  (*»k  ♦  o(£),wA)  ■  (y  ♦  o(j).  Oonaidar  tha  vaetora  a(,)  -  (a^.a^.a^, . ..  ,an>T, 
b(*’  -  (r^.f^ ,...,Fn)*,  and  lat  tha  matrix  I1**  -  l**"1)  ba  dafinad  by 


tB)  |  1  ,  k  -  J  ,  0  <  k,J  <  ■  - 

*kl  1  °hM,J^i  *  *  *  i  «  0  <  k,J  <  ■  - 


.(■) 


aagnituda  of  tha  ia  givan  by 

’•  «i>  °tf t)  -  olirh) 

I 

.<■> 


0  °tT7^  0tTT)  ••• 

o -  : 

o  ~  — o 


.<■> 


IK  I  ia  boondad  aa  n  ♦  •»,  and  wa  can,  for  aaeh  a  >  0,  aalaot  an  a  such  that 


n(m>  I.  4  K.  MOW 


Lat  lb'  l„  <  d,  and 


(I  ♦  *(■,).(■,  -  b(B)  . 
«  <  1.  Hwn 


.<•>  .„<■>  .  i  , 

k-1 


•a<-)  -  «  (1  ♦  l*<")lj"1l*(",l-lb<*ll-  < 


Adding  aoaa  ragularity  aaatnaptiona  on  tha  functiona  to  ba  approx inatad,  wa  can  bound 


(2.19) 


•»  «  -  aup  l»  f I.  , 


If lp  *  sup  |f(z) I 

1  z«r 


Let  f  be  analytic  on  Q,  and  assume  its  Fourier  coefficients  (f(a)  .m*), 

d2* 

w  ■  +(*),  k  -  0,1,2,...,  for*  an  absolutely  convergent  series.  Let  — ?  be  of  bounded 

dw 

variation  on  |m|  ■  1.  Then  is  bounded  wrt  the  norm  (2.19). 

Proof.  He  divide  the  matrix  into  2  parts.  For  an  m  <  n,  let 


A  "  tAkjJ*  Akj  " 


,  k  -  j  or  k  < 


0  ,  else  . 


Let  S  (■  -  A.  For  an  arbitrary  e  >  0,  we  can  select  an  m  so  that  IBI^  <  e  <  1 

for  all  n,  c.f.  the  proof  of  theorem  2.3.  Osing  EA”1  ■  I,  «e  obtain 

<G(n,rt  -  (A  +  I)"1  -  A-1  (I  ♦  *A-1)-1  - 

-  A_1 (I  ♦  t)"1  -  A-1 (I  ♦  £  (-*)k>  -  A_1  +  A*1  l  (-E)k  . 

1  1 


(2.20) 


a  -  A"Vn,f  +  A-1  l  (-l>V*f 

k-1 


sup  I  l  a^p  (s)|  <  sup  |  l  ajwk  +  <  J  l«kl(l  *  «(£))  <  d0u,1 

ser  o  |w|«i  o  o 


for  some  constant  dg  independent  of  n.  From  (2.20), 

lalj  <  (1  ♦  f)"1lA”,l1»lFBfl1  . 
lFnfl.  mas  assumed  bounded  for  all  a. 


COMPUTING  WITH  THE  POLYNOMIALS 


The  proof  of  theorem  2.3  Indicates  a  computational  method  for  determining  the 
coefficients  of  PRf  in  0(n  log(n) )  operations  if  the  Grammian  is  known. 

1)  Determine  n  Fourier  coefficients  b  t"  F  f.  This  requires  0(n  log(n)) 

•  n 

operations. 

n2 

2)  Solve  G  n,a  -  b.  Solutions  of  the  full  system  requires  ~  operations,  but  we 
only  need  to  solve  Aa  »  b,  where  A  is  the  submatrix  of  G*n^  introduced  in  the  proof 
of  theorem  2.3.  A  can  be  choosen  independently  of  n.  Solving  Aa  -  b  requires  0(n) 
operations. 


He  next  turn  to  the  computation  of  the  polynomials.  The  restriction  of  the  mapping 
function  +(w)  to  |w|  «■  1  is  needed,  and  several  numerical  methods  are  available,  see 
Fornberg  [3],  Gutknecht  [5]  or  Reichel  {6}.  The  method  [6]  yields  also  the  capacity  of 
T,  but  not  knowing  the  capacity  only  necessitates  explicit  normalisation 
(Pk(s),(4_1 (s) )k)  *  1,  k  *  0,1, ...,n  -  1.  He  finally  note  that  when  Pnf  has  been 
computed  and  is  to  be  evaluated  at  many  points  it  might  be  advantageous  to  use  a 
representation  which  is  faster  to  evaluate  than  (2.17),  like  a  Newton  polynomial 


representation. 


4.  NUMERICAL  EXAMPLES 

We  consider  two  contours  T,  one  which  is  analytic,  and  ona  for  which  — ■?  has  a 

dw 

jump  discontinuity  on  |w|  m  1.  All  computations  hava  been  carried  out  on  a  UNIVAC  1100  in 
single  precision,  i.e.  with  8  significant  digits.  The  images  of  the  roots  of  unity,  we 
determined  with  approximately  6  significant  digits.  Let  T  be  the  ellipse 
{*  +  iy  i  (j)2  *  y2  “  i}. 

Ex.  El.  j 

S 
10 

Due  to  rounding  errors,  w«  cannot  obtain  a  deviation  much  smaller  than  for  j  »  10. 

Ex.  E2.  Let  f(s)  i<*  (s  *  3)”'. 


1 

sup 

*er 

|(pj+1f)(s)  -  f(s> 

10 

1.3*10 

-6 

20 

3.3*10 

If  the  error  would  decrease  maximally,  see  Gaier  [4],  ch.  1,  it  would  decrease  by  a  factor 
2.56*10  3,  when  j  is  increased  from  10  to  20.  This  is  also  the  case,  when  j  is 
increased  further,  rounding  errors  dominate.  ■ 


Ex.  B3»  Let 

f(s)  «■ 

/*  +  2,  where  we  choose  a  branch  which  has  a  discontinuity  on  the 

negative  real 

axis. 

j 

sup  |(P.  f)(s)  -  f(s)| 

s«r  3 

10 

0.176 

20 

0.121 

40 

0.079 

-10- 


-  «3+1l 


lwl-1 


Ip j+1 (♦<«))  -  w 

4.1*10~ 


2.4*10 


-5 


_ 


I  L  "L 


The  error  decreases  by  a  factor  close  to  — ,  the  expected  rate  of  convergence. 

/ 2 

Zn  the  following  exaaples  r  is  a  sports  ground  shaped  region  obtained  by  placing  a 
unit  square  between  2  unit  disk  halves. 


Ex.  81 


1 

■•*  !?**«<♦<*>> 
|w|-1  1+1 

5 

8.86*10 

10 

2.82*10‘ 

20 

1.24*10 

40 

0.61*10' 

-**1! 

r2 

-2 

i 

-2 

2 


The  error  seeaw  to  decrease  like  o(1/n)f  new. 


Ex.  82  f(s)  i«  (s  ♦  3) 

j 

10 
20 


-1 


ser 


|(P._f)Cs)  -  f(a) | 


J+1 


2.4*10 


-3 


1.0*10 


-5 


Ex.  83  f(s)  i*  /s  ♦  2.  the  same  branch  as  in  ex.  E3 


J 


■ex  |C»4.tf)(s)  -  f(a> | 
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